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Abstract. The INTEGRAL satellite, observing the sky at high energy, has quadru- 
H-5 pled the number of supergiant X-ray Binaries known in the Galaxy and has revealed 

new populations of previously hidden High Mass X-ray Binaries. These observations 
raise new questions about the formation and evolution of these sources. The number 
of detected sources is now high enough to allow us to carry out a statistical analysis of 
PQ the distribution of HMXBs in the Milky Way. We derive the distance of each HMXB 

using a Spectral Energy Distribution fitting procedure, and we examine the correlation 
with the distribution of star forming complexes (SFCs) in the Galaxy. We show that 
f^3 HMXBs are clustered with SFCs, with a typical size of 0.3 kpc and a characteristic 

Mh distance between clusters of 1.7 kpc. 
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1. Introduction 

o 

High Mass X-ray Binaries (HMXBs) are binary systems composed of a compact object, 



CO 



o 



C3 



a neutron star or a black hole candidate, accreting matter from a massive companion 
star: either a main sequence Be star or an evolved supergiant O or B star. Most of these 
sources are observed in the Galactic Plane ( [Bird et aL]|2007| ) as it is expected for such 



young star systems which do not have time to move far from their birthplaces. 
Thanks to the dedicated observations from RXTE and INTEGRAL, around 200 HMXBs 
are currently known in the Milky Way allowing us to focus on their distribution. Using 
RXTE data, |Grimm et al.| ([2002 ) highlighted clear signatures of the spiral structure in 
the spatial distribution of HMXBs. In the same way, Dean et aT] ( |2005| ), |Lutovinov 



eFaLl ( [2005] ), [Bodagh eFeTaL] ( [2007] ) and |Bodaghee et al.| ( |20lTf showed that HMXBs 



observed with INTEGRAL also seem to be associated with the spiral structure of the 
Galaxy. However, the HMXB positions, mostly derived from their X-ray luminosity, 
are not well constrained and highly uncertain due to direct accretion as in HMXB. 
In order to overcome this caveat we present a novel approach allowing us to derive all 
HMXB positions. We build the Spectral Energy Distribution (SED) of each HMXB and 
fit it with a black-body model to compute the distance of each source. Finally, we study 
this distribution and the correlation with Star Forming Complexes (SFCs) observed in 
the Galaxy. 
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2. Deriving the HMXB location within the Galaxy 

In order to compute the HMXB distances, we gathered a set of HMXBs for which at 
least 4 optical and/or NIR magnitudes are known, and for each source its SED was built 
and fit with a blackbody model. This enables us to evaluate the distance of the source 
along with its associated errors. 

2.1. Whole catalogue of HMXBs 

There are currently more than 200 HMXBs detected in the Milky Way. Using an up- 
dated version of the |Liu et aL ( 2006 ) catalogue^ we retrieved the optical, NIR magni- 



tudes, the spectral type and the luminosity class of each source. For each HMXB, four 
magnitudes are required to compute the fitting procedure. Indeed, according to the x 2 
computation, the condition N - n > 2 with N the number of observed magnitudes and 



n the number of parameters left free (in our study, n - 2, cf. section |Z2| ), needs to be 
met. Finally, we selected around 70 sources meeting these conditions. 

2.2. SED fitting procedure 

Our fitting procedure is based on the Levenberg-Marquardt least-square algorithm im- 
plemented in Python. For each HMXB, we build the SED in optical and NIR from a 
maximum of 8 magnitude points (U, B, V, R, I, J, H, Ks) and, as we underlined before, 
a minimum of 4 magnitudes are required by the fitting. This SED is then fit (see figure 
[I} by a black-body model given by the relation: 

A 4 fim -fe 



(: 



where X is the wavelength in yum, Fx the flux density in W.m" 2 .yum _1 , h the Planck 
constant, ks the Botzmann constant, c the speed of light, A A the extinction at the wave- 
length A, R/D the stellar radius over distance ratio and T the temperature of the star. 
From the spectral type and the luminosity class, we derive the radius and the tempera- 



ture of the companion star which dominates the optical and NIR flux (Morton & Adams 



1968[|Panagia|1973]|Searle et al.|20 08). Two parameters are left free: the extinction in 



V band Ay and the ratio R/D whereas the extinction A A is derived from Cardelli et al. 



( |1989| ) at each wavelength assuming Ry = 3.1 for the Milky Way. Knowing the radius 
R of the companion star (see |Morton & Adams|1968[|Panagia|1973| ) , we calculate the 



distance D in kpc. 

The least square function given by the formula 

^i,obs ~ ^i,model 



* 2 =I 



0~i 



(2) 



(with Xi^ b s the observed flux value for the filter /, X^ mo d e u the theoretical flux in 
the i th filter derived from the black-body model and crt, the flux error in the same filter) 
is then minimized by the Levenberg-Marquardt algorithm. 



! see IGR Sources web page maintained by J. Rodriguez & A. Bodaghee (http://irfu.cea.fr/Sap/IGR- 
Sources/|> 
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le-13 Least-square fit to 1E1145. 1-6141 SED 

— Fit 

observed SED 




d = 8.8 kpc 
Av = 5.2 



Figure 1 . Result of the fitting for one source with distance D and extinction in V 
band, Ay. 



2.3. Uncertainties in the computed distances 

The magnitude uncertainties are retrieved from the literature. For the sources for which 
no error is given, we use a systematic error of 0. 1 magnitude. The flux uncertainties 
are then derived from these magnitude errors. Otherwise, we assume the spectral type 
given in the literature as the real spectral type of the companion star. Simulations were 
carried out and enabled us to constrain the spectral type error that appears to be less 
important for supergiant stars than for Zero- Age Main Sequence ones. 

Degeneracy between several parameters values (based on the fitting procedure) 
needs to be taken into account. Indeed, solely relying on a single best fit does not 
capture the full phenomenology associated with SED fitting because extinction Ay and 
distance D are degenerated in this approach. In order to produce the most likely set 
of fits and to determine the dispersion on distance and extinction we carried out 500 
Monte-Carlo simulations for each observed source by varying the photometry within 
the uncertainties. Hence, we generate a random number from a normal distribution (as- 
suming the photometric errors to be Gaussian) contained within the error bars for each 
photometric point so that we build 500 new SEDs derived from the original one. These 
500 new SEDs are subject to the same x 2 statistic computation as the one described 
above. Then, we have an entire set of best fits of parameters (Ay, D) and, we are able 
to plot the distribution in the parameters space, showing the distribution of properties 
derived from these Monte-Carlo simulations and especially showing the dispersion on 
distance and extinction. This dispersion value is taken as the error coming from the 
fitting procedure and the median value of dispersion on distance determination (for all 
the HMXBs under study) is 0.75 kpc. 

There are other sources of uncertainties, particularly the infrared excess of Be 
stars due to their circumstellar envelope generating free-free radiation. According to 
Dougherty et al.| ([l994), this excess should not exceed a mean 0.1 magnitude in J band, 
0.15 magnitude in H band and 0.25 magnitude in K band. However, this value corre- 
sponds to absolute magnitude and therefore can be smaller for sources located far away 
and higher for close ones. To take this effect into account, an estimate of the distance 
and extinction Ay is needed. Since these two values are derived from the fitting pro- 
cedure, we are only able to consider as a distance and absorption estimate the values 
obtained without taking this IR excess into account. This approach is finally equiva- 
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lent to adding a conservative excess of 0.1, 0.15 and 0.25 magnitude to the apparent 
magnitude in the J, H and Ks bands respectively. The results are presented on figure 
01. 
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Figure 2. 1) Green circles represent the initial positions of Be stars whereas blue 
stars represent the source positions taking into account the IR excess. 2) Distribution 
of HMXBs (blue stars) and distribution of SFCs (green circles). Circle radi us repre- 
sent t he different excitation parameter values. The spiral arms model from Russeil 
( 2003| ) is also plotted and the red star (at X = 8 kpc, Y = kpc) represents the Sun 
position at 8 kpc from the Galactic center. 



3. Results: distribution of HMXBs and correlation with Star Forming Com- 
plexes distribution 

We present hereafter the distribution of HMXBs in the Galaxy, obtained with our novel 
approach. The spiral arms model given by |Russeil| ( |2003| ) is also presented (see figure 
[2J2). We want to assess the question whether there is a correlation between this distribu- 
tion of HMXBs and, in a first time, the distribution of Star Forming Complexes (SFCs) 
in the Milky Way (given by |Russeil|2003 ), as it is expected from the short HMXB life- 
time. 



The first approach we adopt is to carry out a Kolmogorov-Smirnov test (KS-test) 
in each axis in order to quantify the fact that the two samples are drawn or not from the 
same probability distribution. We got a value of 0.15 for the X axis, a value of 0.25 for 
the Y axis and a value of 0.31 for the galactic longitude. Even if these values are not 
negligible, suggesting that a correlation between the two samples is likely established, 
part of the information is lost because of the projection on the two axis. To overcome 
this caveat, we propose another method described hereafter. 

We suppose that each HMXB (blue stars on figure [3] 1) is clustered with several 
SFCs (green circles on the same figure). Hence, we can define two characteristic scales: 
a typical cluster size and a typical distance between clusters. Around each HMXB, we 
define several circles with different radius (red circles on figure [3] 1) and we finally 
count the number of HMXBs for which at least one SFC is within the specified radius. 

The number of HMXBs for which at least one SFC is within the specified radius 
versus the circle radius is plotted as the blue curve on figure|3]2. The green curve is that 
expected from chance correlations assuming the HMXBs are evenly distributed across 
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Figure 3. 1) Description of the method used to evaluate the correlation. 2) Result 
of the correlation determination in 2D. 



the sky. The dashed line represents the difference between the two previous curves. This 
difference, being non equal to zero, allows us to state that a strong correlation exists 
between the HMXB and the SFC positions in the Milky Way. Moreover, we compute 
the two characteristic scales exposed before: the typical cluster size of 0.3 kpc and the 
typical distance between clusters of 1.7 kpc, which is larger than the median error on 
HMXB distances. If we take into account the uncertainties in the HMXB positions 



(obtained via the Monte-Carlo simulations) and in the SFC positions (given in |Russeil 
2003, median error of 0.25 kpc), the correlation still exists with the same cluster size 
and the same distance between clusters. These results, while obtained with a different 
method, are consistent with those reached by |Bodaghee et al. ( 201 1[). Fin ally, we test 
our correlation code using a sample of globular clusters ( Bica et al.|2006[ ), principally 
located in the Galactic bulge. Therefore, figure [4] 1 shows the two distributions in the 
Milky Way whereas figure[4j2 shows the result of correlation test. Clearly, as expected, 
no correlation is revealed. 









60 
50 
40 

1 30 

8 

| 

10 



■HMXBs — 


Random — 


Difference 














'■...'■la n -d» * 'c^** * * 

V\x\^ 6>^j4r.. . If. Jj 


|\ 








-15 


-10 -5 5 10 
X (kpc) 


15 


10 
1C 


" : 10 4 


10 D 10 
radius (kpc) 


i 


3 2 



Figure 4. 1) Distribution of HMXBs and globular clusters. 2) Result of the corre- 
lation determination with globular clusters. 
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4. Conclusion 

Evaluation of HMXB distribution is now of major interest in order to study in depth the 
formation of these high energy sources. However, HMXB locations are usually poorly 
constrained and largely dependent on the determination method. We propose here to 
determine the location of a whole set of sources using the same approach: a SED fitting 
of the distance of HMXBs. This method, based on a least-square minimization, enables 
us to reveal a consistent picture of the HMXB distribution, showing them to follow the 
spiral arm structure of the Galaxy. The consideration of uncertainties leads to a small 
error on the source locations and allows us to tackle the study of the correlation with 
SFC distribution. This study shows that HMXBs are clustered with SFCs and enables 
for the first time to quantitatively define the cluster size (0.3 kpc) and the distance be- 
tween clusters (1.7 kpc). The challenge is now to quantify this correlation by taking 
into account the offset between current spiral density wave position and HMXB posi- 
tions expected due to the fact that the matter rotation velocity is different to the spiral 
arm rotation speed. 
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